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Abstract 

We investigate the connection between the well known Sherrington-I 
Kirkpatrick Ising Spin Glass and the corresponding Lattice Gas model 
by analyzing the relation between their thermodynamical functions. 
We present results of replica approach in the Replica Symmetric ap- 
proximation and discuss its stability as a function of temperature and 
external source. Next we examine the effects of first order Replica 
Symmetry Breaking at zero temperature. We finally compare SK re- 
sults with ours and suggest how the latter could be relevant to a 
description of the structural glass transition. 
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1 Introduction 



It is well known [1] that Ising Model is equivalent to Lattice Gas, a system 
denned in terms of occupation variables r taking values and 1. The Lattice 
Gas effective Hamiltonian is formally identical to Ising one[l]. A simple 
change of variables (a = 2r — 1) maps each of the two Hamiltonians into 
the other, provided that Ising external field is related to lattice gas chemical 
potential by h — J = where J is Ising spin-coupling related to lattice 
gas site-coupling $ by J = This results in a simple relation between 
Ising free energy density and Lattice Gas pressure: p = h — | J — / . The 
two systems have therefore the same phase diagram and the same critical 
behaviour (real gases and Ising magnets are in the same universality class). 

For random systems [2] this whole argument breaks down because the re- 
lation between chemical potential and magnetic field involves the quenched 
couplings. As we shall see in the next sections this results in new and unex- 
pected features for the phase diagram of the system. 

For Neural Networks this inequivalence between spin(±l) and occupation 
(0,1) variables had already been pointed out and analyzed, see for example 
[3] and references therein. 

Recently much effort has been devoted to develop a description of the 
structural glass transition[4, 5, 6] within the framework of disordered systems. 
Much of these models were however based on Ising Spin variables instead of 
Lattice Gas ones which should be more appropriate for a condensed matter 
system. For disordered systems the two kinds of variables are not equivalent. 
To have a comparison term it would be useful to analyze the properties of 
a mean field disordered lattice gas model. This could also show us how the 
structural glass transition may be described in the same way as the liquid-gas 
one. 

In section 2 we present our model and in the following one we analyze 
its ground state. In section 4 we discuss its relation with the SK model. 
Section 5 is dedicated to writing down the saddle point equations of the 
replica approach. In section 6 we analyze the low temperature behaviour 
of the Replica Symmetric approximation and show its Phase Diagram. A 
condensed version of section 5 and 6 was originally presented in ref. [7]. 
In section 7 we discuss the stability of the Replica Symmetry and look for 
the AT line of the model. Section 8 contains some preliminary results with 
broken Replica Symmetry and in section 9 we present a comparison of our 
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main results with the SK's ones and draw some tentative conclusions. 



2 The Model 

We consider a system of N sites. An occupation variable r k) defined in 
each site k, can take the values or 1. The Hamiltonian of the system is 
taken to be formally identical to the SK's one[9, 10]. The interaction energy 
between two different (k and /) occupied sites is taken to be <fiki and the 
system is coupled to some external source g. The total effective Hamiltonian 
is therefore: 

N x 

h$[t\ = -g^n- -^2 (f) k iT k Ti . (l) 
k=i 1 i^k 

In magnetic language g would be the external field, while for a Lattice 
Gas g is the sum of the chemical potential, the kinetic contribution, and 
eventually an external force term. 

The infinite-ranged interaction energies {(pki} are taken to be quenched 
independent Gaussian variables with zero mean and variance $ 2 /iV. This 
means that the model should not have, on average, a preference between 
an ordering transition in two sublattices, as in binary alloys, and gas-liquid 
transition. For the low-temperature behaviour we could expect al least two 
regimes. For a strong enough negative g occupied sites should be energetically 
penalized and we expect, for most of the realizations of the noise, the ground 
state energy to be zero and the low temperature behaviour to be similar to 
that of a pure system. On the other hand, for a strong enough positive g 
occupied sites should be energetically favoured and we could expect the low 
temperature behaviour of the system to be much similar to that of SK's one. 
In the following we shall take $ as the unit of energy and set $ = 1, so we 
have: 



p (4>ik) = [tt) ex P 



,2tt7 "' r I 2 TlK \ ^ 
Each (pki is taken to be equally distributed and therefore each site inter- 
acts with each other. We also set the Boltzmann constant equal to 1, as a 
consequence T, H, g and <fi are all dimensionless. 

For a given realization of the 0's the partition function is: 

ZM9) = H^ mT] - (3) 

{r} 
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For a Lattice Gas this would be the Grand Canonical partition func- 
tion, because the total number of occupied sites J2k T k i s n °t constrained to 
a specific value. Strictly speaking we should call "pressure" the thermody- 
namic potential (InZ), while the grandcanonical independent variable should 
be the chemical potential. To have a more immediate comparison with SK 
model, and also with other models based on spin variables, we shall however 
abusively refer to the thermodynamical variables as if they were the canon- 
ical ones. In the following we shall therefore call "free energy density" the 
thermodynamic potential and refer to g as the "external field" . 

We are interested, as usual when dealing with quenched disorder, to eval- 
uate the averaged free energy density: 

f(P;g) = -^ J P[0]in^#, (4) 

this will be done in the following using the replica approach[2, 8, 9]. 

3 Ground State Properties 

We now discuss in more detail the ground state picture which was conjectured 
in the previous section. First of all we observe that the energy of the con- 
figuration without occupied sites is zero for all realizations of the quenched 
couplings. We can therefore conclude that the ground state energy is always 
non-negative. Secondly, we can easily see that the configuration with all oc- 
cupied sites has an average energy equal to —gN, which is exactly the same 
value obtained for the all-spins-up configuration of SK model. Combining 
these two results we obtain, for the (quenched averaged) ground state energy 
density, the following upper bound: 

u < -g6(g) 

A lower bound can be obtained from the knowledge of the eigenvalues for 
a large Gaussian random matrix. Let us call uo the maximum eigenvalue of 
the matrix [4>ik\- For every configuration we can write 

ik k 
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In our case we have {rk) 2 = Tfe, and u = 2 with probability one in the 
thermodynamic limit. Putting all together we obtain, for the Hamiltonian 
(1), the following bound: 

-(i + s)E^<^M • 

k 

Taking the minimum of both sides we finally get: 

-(l + g)e(g + l)<u. 

We thus see that for g < —1 the upper and lower bounds saturate and 
u — with probability one in the thermodynamic limit. The ground state 
configuration is that with no occupied sites. Each configuration with a finite 
number of occupied sites has the same average energy density in the thermo- 
dynamic limit. We therefore conclude that the low temperature behaviour 
of the system should be very similar to that of a pure one. 

On the other hand, for positive g we have u < —g, bound which is 
saturated in the limit g — > oo. This implies that our system has a (first 
order?) zero-temperature transition at some go between —1 and 0. We also 
speculate that for a large positive g the low temperature behaviour of the 
system should approximate that of SK model. We note that in SK case 
repeating the previous arguments we would obtain —1 — \h\ < u < —\h\. 

The lower limit for go can be improved using Derrida's argument [11]. Let 
us define T,(g,N) the average number of configurations with negative energy 
for a system with N sites: 

J T 

For a negative g we obtain, in the thermodynamic limit, 
a= lim ±-\nZ(g : N)=\n2-g 2 . 

If g < — Vln 2 we see that a < and therefore N) — > in the thermody- 
namic limit. Following Derridafll] we can conclude that for g < — Vln2 = 
—0.8326 there are almost surely no configurations of negative energy in the 
thermodynamic limit. Since there always is at least a configuration of zero 
energy, this must be ground state energy. We thus obtain the lower bound 
go > — Vln2 = —0.8326 for the transition field at zero temperature. 
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4 Relation with SK model 



Before proceeding in our analysis we want to discuss the connections be- 
tween the SK model[9] and the previously introduced one. We assume, as 
for the corresponding homogeneous models[l], that site variables are related 
by: 

<7 = 2r-l <=► r= l -{a + l). (5) 
Substituting (5) into the Hamiltonian (1) we get: 
a+ 1" 



-H, 
where 



\ E {? + ^) + E Q# + ^) ^ + 1 E <fe<^ , (6) 



= E 0« • 



We stress that (6) connects a homogeneous (site-independent) field system 
to a local- field one. In fact, defining the SK effective Hamiltonian for a 
site-dependent magnetic field: 



k Z k^l 



we can write: 



Ha 



(7+1 



1 



and the local magnetic field is correlated to the couplings by: 

h k [(j>] = 2g + $ k . 



(7) 



(8) 



(9) 



The change of variables (5) maps our Hamiltonian (1) into an SK-like one, 
but correlation given by (9) is enough to destroy their equivalence. In fact, 
summing over configuration, we get for the partition functions: 



Z^ifog) = exp 



(10) 



6 



and thus the relation between the free energy densities of the two models 
leads to: 

c(SK) ( P. 



mg) = -\g + \ffi, 



(ii) 



Our system is therefore equivalent to an SK-like one in which the magnetic 
field is a local random variable correlated with the couplings. SK's averaged 
free energy is not directly related to ours. Relation (5) is thus not useful in 
investigating thermodynamics of our model that is not reducible to something 
known, we have to face it by itself. 



5 The Replica Symmetric Solution 

Let us now proceed in applying replica formalism [2] to our system; we have 
to calculate the averaged n-th power of the partition function: 



z n = {z^y = J {z^p[<t>]d(t> . 

For integer n we get, after performing Gaussian integration: 



z n = E ex P 



e/^e^+^e(e« 

k a, ZJV k<l \ a 



using the identity: 



(12) 



2E E« =E E 



r k r k 



k<l \ a 



a,b \ k 



EE^ fe 

k a,b 



(13) 



which follows from the relation r 2 = r, we can reorder the exponent and 
obtain: 



Z n = E ex P 

{t} 







k \ a a ,b 



_ab 



IJexp 

a,b 



4 fP 



Nf3 2 



-rll T k 



ah 



Using Gaussian identities we rewrite (14) as: 



(14) 



Z n = E eX P 

{r} 



P 



/^EUE^-^E^ 

k \ a a,b 



X 
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a,b 



'N(3 2 
An 



exp ^-j(3 2 Q 2 ab + \j5 2 £ Qa b T a k r£j dQ ab 



reordering the exponentials and defining: 



a,b 



.M 



we finally get: 



s„(/3;<7) 



'N(3 

47T 



2\ — 



The averaged free energy density is given by: 



(15) 



(16) 



(17) 



(18) 



f(P; 9) = J im lim 



TV^oo rwO /3niV 



In Z n (/3;g) 



(19) 



In the thermodynamic limit the integral can be estimated by maximizing the 
integrand, and this yields: 



fn(P;g) = Jim 



1 



lnZ n (P;g) 



1 



n^oo pnN "^'^ j3n 



M{A[Q]} . 



(20) 



The extremum is determined from the saddle point equation: 

OA & p 2 J2 T r a r b e~ H ^ 
dQ7 b ~ Y Qab ~ T E.e-^M " ' 



(21) 



that may be rewritten as Q ab = (T a T b )Q. 

We first consider saddle points that are symmetric under the Replica 
Group [2]. Setting Q ab = q + b5 ab we can write J2abQabT a r b = g(Ea r<1 ) 2 + 
bJ2a T a and J2 a bQlb = n (<? + b) 2 + n(n — l)q 2 , substitution in (17) and ex- 
traction of the n — * limit then yields: 



1 1 poo r 

/ = -8b(2q + b) - (2tt/3 2 )" / In 1 + S a+ *>® 
4 J— 00 L 



e 2 



(22) 
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In equation (22) we set a = g + \(3b, and the matrix elements satisfy the 
coupled equations: 



_1 -±z 2 1 

e 2 



p = q + h= (2tt)-3 /" [l + e /3(a+z v / 9)" 

g = (27r)-3 y [l + e^+^j ~ 2 e -^ 2 rfz . (23) 

As can be seen following the line of [9, 10] the physical significance of p 
and q is: 

P = F) Q = W, (24) 

where, following the notations of [2], a bar denotes the average over quenched 
disorder. 

For (3 = we have p — \ and q — |, as we expect from their physical 
significance. In the high temperature regime we can solve (23) by expansion 
in powers of f3, and this yields: 



6 Low temperature results 

At zero temperature we can perform a detailed analytic study of saddle 
point equations. In this limit eq. (23) leads to: 

q = (2tt)-3 r /v V^ 2 dz = erf f— ) (25) 

7o = lim P(p -q) = (27rg)~^ e ~% , 

the equation for p is the same as that for q indeed for T — > we have 
p = q + 7T and thus a = g + |7o- 

For g < there always is a solution of (25) with q = 0, 70 = and a = g. 
To look for solutions with g 7^ let us define x = erF x (g) = a/y/q, then we 
get: 

p = q = eri(x) 
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Figure 1: External field g versus parameter x (x = erf (q)) as results from 
the last of (26). 
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g 



-0.8 -0.6 -0.4 -0.2 0.2 

Figure 2: Free energy density versus external field at zero temperature in 
the replica symmetric solution. The transition point is where the two lower 
branches cross each-other. The higher branch is an unstable solution. 



7o = [2tt erf(x)] 2 e ^ 
a = x\J erf(x) 

9 = a - = %\j erf (a;) - - [2vr erf(ac)]~ s e~^ 2 , (26) 

we can therefore express all relevant quantities as functions of parameter x. 
From the last of (26) 's (see fig. 1) we can see that each value of g in the 
range —0.70242 < g < corresponds to two values of x both lower than 
0.30859. Summarizing, if g < —0.70242 or g > we have a single solution of 
(25) for each value of g, but if —0.70242 < g < we have three solutions and 
in order to pick-up the physical one we have to impose the continuity of the 
free energy density as a function of external field, see fig. 2. We thus find a 
first-order phase transition in the point go = —0.63633 where the two lower 
branches of / cross each-other, the transition point being determined by the 
condition 70 (go) = —go- We show in fig. 3 and 4 q (= p) and 70 as functions 
of external field. As expected both are discontinuous on the transition point. 
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-1 -0.5 0.5 1 1.5 2 2.5 3 

Figure 3: Replica symmetric order parameter q (= p = (r) ) at T = as a 
function of external field. 

Next we look to thermodynamical functions, the internal energy and entropy 
densities are given by: 

u = -(g + lo)q 
s = ■ (27) 

The entropy (27) is negative in the range g > go where 70 is different 
from zero, so we should expect the Replica Symmetry to be broken in this 
region. We stress that, differently from SK case, we have a region in which 
the Replica Symmetric Solution remains physical down to zero temperature. 
The maximum absolute value of the zero temperature entropy is at g = g$ 
where it takes the value 0.101 and it strongly decreases for higher values of 
g (e.g. at g = 1 s = -0.011). 

For /3 ^> 1 the solution of (23) can behave in two different ways. In the 
range g < g , where q (= q(T; g)\ T=0 ) and 70(g) are identically zero, we 
find that all their temperature derivatives vanish for T — > 0, q vanishes as 
e /3g < 0) anc [ ry (= p _ ) as fj^Q _ Otherwise if g > g , they depend 
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Figure 4: Behaviour of 7 (= lim T ^ (3{{r 2 ) — (r) 2 ) ), in the replica symmet- 
ric solution, clS db function of external field. 
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0.4 0.8 1.2 1.6 2 



Figure 5: Replica symmetric order parameter p as function of temperature 
for several values of external field. From top to bottom the values of g are 
0.5, 0, —0.6, —0.68, —0.7 and —0.75. The embedded picture represents a 
magnification of the region around the transition. 



linearly on T, indeed we get: 



= go7o(go + |7bc*o) t 

° (go + |7o«o) 2 + i7o(<7o - (A) 



ho(9b-«o) 



7 7 ° + (<7o + !7oao) 2 + i7o(<?o - al) T 
p = q + 1 T . 

We have also numerically solved (by iteration) equations (23) for several 
values of T and g, results are plotted in fig. 5 and 6. We find a line of 
first-order phase transitions in the phase diagram of fig. 7. Such line ends in 
a second-order transition point for T ~ 0.22 and g ~ —0.7, as attested from 
fig. 8 where the linear response function x — dp/dg (the "susceptibility") is 
plotted versus temperature for g = —0.7. 
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Figure 6: Replica symmetric order parameter q as function of temperature 
for several values of external field. From top to bottom the values of g are 
0.5, 0, —0.6, —0.68, —0.7 and —0.75. The embedded picture represents a 
magnification of the region around the transition. 
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Figure 7: Phase diagram for the replica symmetric solution on the plane 
(g; T). A line of first order phase transitions ends with a second order tran- 
sition point. 
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Figure 8: Linear Response Function \ — dp/dg ("susceptibility") as a func- 
tion of temperature for g = —0.7. 
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7 Study of stability 



As we pointed out in the previous section the Replica Symmetric ansatz 
for the saddle point equations (21) can lead to unphysical results at low 
temperature. This should be considered gnal[2] of Replica Symmetry 

Breaking (RSB). To check this conjecture we have to calculate the eigenvalues 
of the Hessian of (17), this will be done in this section following the line of 
ref. [12] but with two main differences. The first refers to the fact that we 
have = but the diagonal part of the matrix Q ab is different from zero 
and plays the same role as the magnetization vector of ref. [9, 10, 12] where 
J = J 7^ 0. The second difference is that we never impose a priori Q ab to be 
a symmetric matrix. As a matter of fact the only term in A[Q] depending 
on anti-symmetric part of Q ab is Tr Q 2 , its contribution could be integrated 
out and has no physical significance but allowing for its presence permits a 
more compact and elegant notation. We can indeed treat the diagonal and 
the off-diagonal part of Q ab in the same way. However we shall see that 
the saddle point is (as it should be) always a symmetric matrix and stable 
against anti-symmetric perturbation. 
The Hessian of A[Q] is: 



TTCd 

H ab = 



f) 2 A 1 

dQ^dQ7 d ~ Y dacdbd + ab ' 



where we have set: 

and the elements of V have the following symmetries: 



-ped -pea -pac -p 

ab ba ab cd 



^cd 



^dc 



ab 



Remembering that r 2 = r, and substituting the Replica Symmetric ansatz 
for the saddle point, we have the following different matrix elements for T: 



1 



1 



r c 



\(3 4 [(r a ) Q (r c ) Q -(r a r c ) Q ]= 1 -(3 i (p 2 - q ) 
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Kl = -/ [(r a r%(r a ) Q - <rV%] = \^ (qp - q) 



1 



rS = ~f [(r a r b ) Q (r c ) Q - (AV) ( 



r/3 4 (?P " A*) 



ab — 



t r 



= \p A [<tV> q <tV% - (rW%] = (g 2 - p 3 ) 

= X ~/ [(A%(r c r% - <rVrV%] = (g 2 - p 4 ) , 

where explicitly written replica labels are different. In the n — > limit, the 
required expectation values are: 

(T a )Q* =QL = P 

(T a T%* =Q* ab = q 



(rVr%. = p 3 = (2tt)- 



1 + e 



(tVtV%* = p 4 = (2tt)- 



]_ _|_ g-ZS^+^v^) 



e 2 2 ^ 



1-4 12 



e 2 Z (iz. 



The eigenvalues equation for the Hessian if^ : 

H abVcd = ^Vab , 

c,d 



(28) 



has, for general n, four classes of eigenvectors with no more than six distinct 
eigenvalues. Anti-symmetric eigenvectors (i] ba = —i] ab ) give the eigenvalue 
— \l3 2 which (as we should expect) is always positive. Eigenvectors that 
are symmetric matrices, and invariant under the Replica Group, give two 
eigenvalues that, in the n — * limit, are: 

Ai, 2 = i/3 2 +^ 4 (8p 3 -6p 4 -g-p)±^ 4 v /(8p3 - 6p 4 - 3g + p) 2 - 16(p 3 - q) 2 - 

(29) 

Symmetric eigenvectors (r] ba = r) ab ) that are invariant under interchange of all 
but one of the replicas give, for general n, two more eigenvalues A 3 4 that for 
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n — > reduce to the previous ones. There are finally symmetric eigenvectors 
that are invariant under interchange of all but two replicas, for n — > these 
give rise to the eigenvalue: 



(30) 



In the high temperature regime the eigenvalues Ai and A 2 are found to be, 
in spite of the hermiticity of the Hessian matrix, complex conjugate. This 
should not be too surprising because we are working in a space where the 
norm is not positive definite (lim n _>o - Tr Q 2 = p 2 — q 2 ). However, as could 
be seen considering the gaussian approximation to the integral in (18), we 
believe that the stability of the solution is determined by their (common) 
real part which is always positive. At lower temperature these eigenvalues 
are real and never negative, although they are found to vanish linearly on 
the second-order transition point. This means that this saddle-point is stable 
against Replica- Symmetric perturbation. 

The eigenvalues Ao A3 and A4 are not relevant and thus the stability of this 
solution against Replica-Symmetry-Breaking perturbations is determined by 
the eigenvalue A 5 . This eigenvalue is always real, to study its sign it is useful 
to define A 5 = 2T 3 A 5 , so we get: 



The low temperature limit then follows in a straightforward way and we find: 



It is therefore apparent that, if g > go, A5 must become negative at low 
enough temperatures. The Replica- Symmetric saddle point is unstable and 
the Replica Symmetry is spontaneously broken. We plot in figure 9 the be- 
haviour of A5 versus temperature for several values of g. The line of insta- 
bility obtained by numerical evaluation is shown in figure 10, we stress the 
presence of a region in which the Replica Symmetry remains exact down to 
zero temperature. 



(a + 2Tu) 2 




(31) 



lim A 5 = — (2irq) 2 e 2 q 
t^o 6 



1 
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0.1 0.15 0.2 0.25 0.3 0.35 0.4 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

Figure 9: Behaviour of A5 (= 2T 3 As) as a function of temperature for several 
values of external field. From bottom to top the values of g are —0.6, —0.68, 
—0.7 and —0.75. The embedded picture represents a magnification of the 
region around the transition. 
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8 Zero temperature results with first order 
Replica Symmetry Breaking 

Once stated that the Replica Symmetry can be broken, we looked for saddle 
points that are not invariant under the Replica Group. As a first step we 
analyzed the zero temperature behaviour of the solution given by the ansatz 
proposed in ref. [13, 14]. Let us divide the n replicas in v = n/m groups, 
each of m replicas, and take the matrix elements of Q as follows: 

Qab = p = q + d + b if a = b 

Qab = <li = <l + d if a 7^ b but I(a/m) — I(b/m) 

Q ab = q = q if I(a/m) + I(b/m) , (32) 

where as usual [13] I(x) = mm{n e N : n > x}. With this position we have 
Tr Q 2 = n [p 2 + (m — l)qf + (n — m)q' 2 ] and the effective Hamiltonian (16) 
reads: 



H Q [r] = -(3a^r a - -J\, fe^V " ^E ( E 

a Z \ a / z i \jfc( a )=i / 



(33) 



where fc(a) = I(a/m) and as previously we have set a = g + 

Using standard properties of Gaussian integrals, and setting g a (x) 

_1 

(27ro")~2e~2^, we find in the n — > limit: 



where: 



Imy^i d) 



9d{y)dy . 



(35) 



Substituting the ansatz (32) in the saddle point equations (21) we find, in 
the n — > limit, the following equations for p, qi and qo: 



P 



J g qo (z - «) 

J g qo (z-a) 



I 



1 + e 



m—1 



P(z+y) 



9d(y)dy 



I m (z; d) 



dz 



\ _|_ e P(z+y) 



m—2 



9d{y)dy 



I m (z; d) 



dz 



(36) 
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/ fl + e^ +!/ ) 



m—l 



9d{y)dy 



I m (z; d) 



dz 



To decide if the correct saddle point has to be a maximum or a minimum 
of eq. (34) we note that lim n ^ ^ Tr Q 2 = p 2 — (1 — m)q\ — mq 2 so, for 
< m < 1, we should expect the saddle point to be a minimum with respect 
to p (diagonal parameter) and a maximum with respect to qi and go (off- 
diagonal parameters). Because m itself is a parameter for the off-diagonal 
part of Q a b we speculate that / should be maximized with respect to it [2]. 

In order to work out the T — > limit of eq. (36) let us consider the 
internal energy density: 

u = -gp - ^/3b(p + gi) - ^pmd(q + q ± ) . 

To keep u finite as (5 — > oo we set b = 7T and m = pT, and we expect that 
7 and p remain finite as T — > 0, indeed in this limit eq. (36) leads to: 



z)g qo (a + z) 



7 ' erf(^)+ e -^-^erf(^) 



dz 



q + d = J 



1 



erf(^)+e-^-^) erf(^) 
e-**-\rt erf(^) 



g qo (a + z)dz (37) 



erf(^)+e-^-^) erf(^)J 
In the same limit the free energy density (34) then becomes: 



where 



/ = ^ [2 7 (g + d) + pd(2q + d)} - - J g q (a + z) In [I m (z; d)} dz 

pd — z\ 



lim I m (z; d) = erf + e -"(*-§/«0 erf 



(38) 



(39) 



Equations (37) have been numerically solved by iteration with p held fixed, 
next the resulting free energy (38) has been maximized with respect to p 
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Figure 11: Free energy density versus external field at zero temperature with 
first order replica symmetry breaking. The first order transition point is at 
g = -0.625. 
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Figure 12: Behaviour of /x (= lim^o f3m) as a function of external field. 
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Figure 13: Solution with first order replica symmetry breaking: Off-diagonal 
order parameter q (= q ) at T = versus external field. 
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Figure 14: Solution with first order replica symmetry breaking: Symmetry 
breaking parameter d (= q% — q ) as a function of external field at zero 
temperature. 
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w - w -> — I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-1 -0.6 -0.2 0.2 0.6 1 

Figure 15: Solution with first order replica symmetry breaking: Behaviour 
of 7 (= liniT_ i .o /3(p — qi)) as a function of external field. 

using a standard IMSL routine. Such numerical solution of the saddle point 
problem gives the results plotted in figures 11-15 as functions of g. We note 
that the transition point is shifted from g Q = —0.63633 to go ~ —0.6250 and 
the discontinuity of q (0.27) is split in 0.28 for p and 0.22 for q ; we also note 
that the values of 7 are about a factor four smaller than those obtained in 
the Replica Symmetric approximation. 

9 Discussion conclusions and outlook 

In this paper we have investigated the thermodynamical properties of a 
lattice-gas model with infinite-ranged random interactions. We have dis- 
cussed its relation with SK model and showed that the averaged thermody- 
namical functions of the two systems are not directly related to each other. 

In SK case the zero-field Hamiltonian has a global Z2 symmetry and all 
thermodynamical functions are either even or odd in the external field. SK 
model can have a first-order phase transition, as a function of magnetic field, 
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only if a strong enough ferromagnetic part (i.e. a non-zero mean) is added to 
the random coupling, such a transition is indeed related to the breaking of the 
global Z 2 symmetry induced (as in the homogeneous case) by a ferromagnetic 
coupling. Moreover, at each value of external field, SK's Replica Symmetric 
Solution always becomes unstable for low enough temperature. 

It is apparent how our picture is different from the usual one. We have 
considered the case of purely (zero-mean) random interactions and our Hamil- 
tonian has no Z 2 symmetry. Nevertheless the Replica Symmetric Solution of 
our system exhibits a line of first-order phase transition points ending with 
a second-order transition point. This feature is robust to Replica Symmetry 
Breaking, indeed the second-order transition point lies outside the unstable 
region and we only expect that the transition line should be deflected on the 
instability boundary. Two phases co-exist along the transition line and, as 
we have stressed, in one of them the Replica Symmetry is exact down to zero 
temperature. We think that this is a consequence of the fact that in this 
phase the system should behave like an homogeneous one. 

To have a physical interpretation of the obtained phase diagram it is use- 
ful to return to a more conventional framework for lattice gases. The effective 
Hamiltonian (1) describes a system of particles, in a discrete space, interact- 
ing with a two body potential <pki- In this context g is closely related to the 
chemical potential and the thermodynamic potential (4) is the negative of 
the pressure. The parameters p and q are respectively the quenched averages 
of particle density and of its square. When a first-order transition occurs 
there are two coexisting phases, the low-density one is to be interpreted as 
the gas phase, while the other is the liquid one. In our model the gas phase 
is that which is replica-stable at all temperatures. On the other hand, at 
sufficiently low temperatures, the Replica Symmetry breaks inside the liquid 
phase giving up to a glassy state. A part of the line of first-order transitions 
may therefore be interpreted as glass in equilibrium with its vapour. 

We thus have a simple and soluble mean-field model accounting (in a 
conventional and so far well understood way) for both a liquid-gas transition 
and a glassy regime. This could be regarded as an intriguing paradigm for 
capturing the structural glass transition in the framework of replica theory. 

There is a flaw in this. We have analyzed a model with quenched (ran- 
dom) disorder included by hand. Structural glasses instead do not necessarily 
have random interactions in their hamiltonians. In several recent publica- 
tions (see e.g. [4, 5, 6]) the main effort was first to present a system, without 
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random interactions, which behave in a glassy way, and next try to construct 
a disordered model which mimics the starting one. Here we have to go in the 
opposite direction. We have a good candidate of random system, the way of 
proposing a corresponding deterministic one will be presented elsewhere[15]. 
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